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Abstract 

We analyze the asymptotic behaviour of the heat kernel defined by a 
stochastically perturbed geodesic flow on the cotangent bundle of a Rieman- 
nian manifold for small time and small diffusion parameter. This extends 
WKB-type methods to a particular case of a degenerate Hamiltonian. We 
derive uniform bounds for the solution of the degenerate Hamiltonian bound- 
ary value problem for small time. From this equivalence of solutions of the 
Hamiltonian equations and the corresponding Hamilton Jacobi equation fol- 
lows. The results are exploited to derive two sided estimates and multiplicative 
asymptotics for the heat kernel and the trace. 



1 Introduction 

The problem of finding Gaussian bounds for the heat kernel and of deriving asymp- 
totic expansions for non degenerate diffusions has been studied by various researchers. 
Expansions for the stochastic flow and stochastic variational calculus have been 
exploited by |Mol75| . |Dav88| . |ABdMBB95] . |KS87| . |BA89j . and |Lea94| . For 
hypoelliptic operators satisfying a uniform weak Hormander hypothesis uniform in 
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space estimates for the heat kernel have been given, e.g., by Leandre. For Gaus- 
sian hypoelliptic generators D. Manankiandrianana |CME81| and M. Chaleyat, L. 
Elie |Man79| have classified the structure of the corresponding processes, derived 
the associate action functionals, and the associated subriemannian metrics. 
In |Mas85| Maslov observed that for a class of pseudodifferential equations (called 
tunnel equations) on M n which includes the example of elliptic second order equa- 
tions asymptotic bounds for the heat kernel can be established for small values of a 
parameter h (Planck's constant). Results using WKB-type methods may be found 
e.g. in |Mas72| . |Mas85| . |MF81| . |KM97| . |KolOO| . In this paper we modify the 
WKB-type method by introducing a change of variable which adjust the different 
time scales of the variables. We derive bounds for the kernel of degenerate diffusions 
associated with the geodesic flow on a compact <i-dimensional Riemannian manifold 
M with metric g which is subject to a stochastic perturbation in the tangent space 
variable. In local coordinates the diffuion on the cotangent bundle T*M is given by. 

= : hT x B. ( x,y, -h^-, -h-^- ) u 
\ ax ay J 

with initial data u(0, x, y, xo, y°) = 5(x — xo)5(y — y°) for x and y in a neighbourhood 
of the point (xo,y°) G T*M, h > is a (small) parameter gij(x) is the expression 
of the metric in local coordinates with inverse G^ix). The summation convention 
over indices is used. The pseudodifferential operator H(x,y,—-^,—j^) defines a 
Hamiltonian function H(x,y,p,q) the arguments of which can be divided into the 
variables x, y G M. d and their respective conjugate variables i.e. the momenta p and 
q. In local coordinates the Hamiltonian takes the form 

H(x,y,q,p) = ^ (g(x)q,q) - (G(x)y,p) + ~ (-^-(G(x)y,y),q\ . (2) 

It is possible to write a stochastic differential equation associated with the operator 
C One way is to introduce an imbedding r of the manifold M into a large Euclidian 
space WL N (such an imbedding always exists). In local coordinates the corresponding 
stochastic system will have the form 

dH , 

axi = — — at 

9H , d , , , 

Vi = ~dx~ dx~^ 1: ' ^ 

where (x.y) G T*M, u is standard Brownian Motion in M. d and (■,■) denotes the 
standard scalar product in WL N . The system describes a Newtonian particle corre- 
sponding to a Hamiltonian system perturbed by a Gaussian White Noise force u w u 
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with the friction \^ l y k yi see e.g. |Nel67| . |AHZ92| . |MSnnj for the flat case or |.Toe78| 
for a covariant definition. In the sequel neglecting friction we shall restrict ourselves 
to the example of the geodesic flow on the cotangent bundle T*M of a Riemannian 
manifold M. In this case the Hamiltonian corresponding to the deterministic part of 
([3]) is of the form f(x) = ^(G(x)y,y) where G(x) = is smooth by definition, 

see also ([2). 

The content of the paper is organized in five sections. The second outlines the general 
structure of an invariant diffusion. In the subsequent one we study the underlying 
Hamiltonian dynamics of the stochastic geodesic flow. This results in the uniform 
estimate (fT9l) and in the existence of a local diffeomorphic flow on the cotangent 
bundle. In Section 4 we explain the WKB-method for imaginary time and the change 
of variables which compensates for the degeneracy and allows for a closed system of 
algebraic equations for the coefficients of an asymptotic expansion for the solution of 
(PQ) in powers of h. The remaining part of this section is dedicated to the existence of 
a multiplicative error. The paper is concluded by discussing an asymptotic expansion 
for the trace formula for the generator corresponding to a stochastic geodesic flow 
(PQ). The coefficients of the expansion are geometric invariants of the manifold. 



2 Geometric background for the system under con- 
sideration 

Let us introduce a diffusion process on the cotangent bundle of a compact Rieman- 
nian manifold M with metric g which in local coordinates has the form 



d h? ■ ■ du 

h-u = -g ij{ x)—u + h(a^x) + ^(x)y 3 ) — 

+ h{bi(x) + P{(x) yj + ^i l (x)y jyi )p- - V(x, y)u (4) 

where h is a strictly positive parameter, the d x c/-matrix g(x) = (gij(x)) is posi- 
tive definite, a, j3, 7 are d x <i-matrices, a being non-degenerate, a, b are uniformly 
bounded in x G M d and V is a polynomial of degree 4 in y G M d with coefficients de- 
pending on x G M. d bounded from below. This diffusion defines a regular degenerate 
diffusion of rank one in the sense of [KolOOj . In this book there are given necessary 
and sufficient conditions for this diffusion to be invariant. 

Let us show that the pseudodifferential operator in (JXJ) is an invariant degenerate 
diffusion. In order to study invariance under change of coordinates in the manifold, 
we recall that changing the variable x — > x for the diffusion (pQ) induces a change of 
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the momentum variables according to the rule y = This implies e.g. 

djjk dxi du du dx l du dye . . 

+ l^^J ( 5 ) 



dyj dx h ' dx l dx e dx l dye dx 
du du dye d 2 u d 2 u dye dijk 
dyi dye dy { ' dyidyj dy k dy e dyj dy { 



(6) 



for i, j, k = 1, . . . , d. An immediate consequence of ((5|) and is that the terms 
of different orders do not mix and hence may be considered separately. Then the 
invariance of equation ((H) follows from direct calculation using the relations j5j) and 
([6]) and the fact that the metric tensor transforms according to gij(x) = gij(x)^^. 
The calculations in this paper will be carried out in normal coordinates around Xq. 
We are using the definition given in [CFKS87] for which x = and for which the 
matrix associated with the Riemannian metric satisfies detg(a;) = 1, which can be 
achieved by a change of variables and the expansion 



9^) = ^ + ^x k x l + 0(\x\ 3 ). (7) 

These conditions imply 

n 

$>£ = for l<k,l<d 

i=l 

and the following representation for the Gausssian (scalar) curvature 

i.k 

For the corresponding inverse matrix G(x) = g^{x) and its derivative we get the 
following Taylor approximations 

G ij (x) = S{- l -glx k x l + 0{\x\') and T$(x) = -g%x e + 0(|x| 3 ). (8) 



2" 

Here T^(x) := ^pr(x) are the coefficients of the Levy Civita connection in local 



coordinates. For symmetry reasons we also introduce the notation juj = which 
will be used in the next section. 



3 Hamiltonian Dynamics 

The Hamiltonian system corresponding to ([I]) is given by 

377 

* = -T- ( - J ^ (9) 
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" " Gv(x) Pj - ^(x) Vj , q = G*{x)pj - r^x) yj q k 



dyi \dx 

dH ( d \ \ f d d \ 1 / d 

with initial condition x(0) = Xo, y(0) = y°, g(0) = go, p(0) = p°- The corresponding 
flow will be denoted by (X, Y, Q, P). 

Hypothesis 3.1. Let there exist constants k,to and Co such that for all c G (0, Co] 

and t G (0,t ] the solution to (Tj|) exists on [0,t] whenever the initial data y°,qo,P° 
are satisfying \y°\ < |, |go| < |p°| < % and xq may be chosen freely. 



Proposition 3.1. Under the hupothesis \3.1\ the components of the flow (X, Y, P, Q) 
started at (xo,y°,p°,qo) satisfy 

\x(t) -x \< kt(l + j), \y(t) - y°\ < kt(l + °-) 

w) - qo \ < kt (i + \ P {t) - P °\< M{i + ° [ 



Proof. The proof uses Taylor expansion and is a direct consequence of the fact that 
the matrix g is a bounded function of the variable x (c.f. also Proposition 2.3.1 in 
|Kol00| V □ 

If we associate with the variables x, y, q,p the degree 0, 1, 2, 3, respectively each right 
hand side of (JH) is of specific degree varying between 1 for x and 4 for p. The degree 
is related to the power of | which is crucial for the iteration procedure giving the 
expansions intrinsic to the WKB-method. 

In the subsequent proposition we show that the solution (x,y,q,p) of ([9]) starting 
at (xo, y°, qo,P°) can be represented in terms of the solution (x,y,p,q) solving 
with initial value (x , y°, 0, 0). For reason of space we shall skip the arguments of 
the functions g(x), G(x) and of their derivatives in future. 



Proposition 3.2. Under the assumptions of Proposition [Q we have: 
x = x-lt 2 qo-^l P +{n>y°)q }+6 4 

y = y + t q + h 2 [p +(Qy )q ] + j5 4 

q = q + t[G{xo)p o + {Ty o )qo] ^ + ^ [iAy oy + ^ 

p = p + tS . 
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where the error term 5 is of the order 0(t + c). Moreover, x = Xq — y^t + 0(t 2 ), 
y = y° + 0(t), p = 0(t), q = 0(t), and 

I J- VoJi - h Vj - ^faTVj ' \ y,y ) ik ~ ikViVj ~ lhF VeVj 

( Q Vo)u = !l,,y J , ! "<lji ~ lujG 3m ] (x )y° m , 

(n'y°) ik = T^g ok y«- llkm G^- gi] TU 

(Ay f = -T ik G mi y° e + G im T^y° e - G m/ T^°, 



1 1 

/A'/\' _ /"Wi,., „ rim 



Proof. For the proof we expand the components p — p, q — q, y — y,x — xin Taylor 
series with respect to time. For the time derivatives of the expansion we insert the 
explicit expressions given by the system of Hamilton equations ([9]). Due to linearity 
of differentiation it suffices to study the derivatives of y(t) and insert go = and p° = 
to retrieve the derivatives of y(t). We begin with the function Api(t) = Pi{t) —pi(t) 
because it is of highest degree and hence determines the degree to be considered for 
all other components. We have Api(t) = p° + J Q Apj(r) dr. We directly find: 

AR(t) = P° + jf if yePk - \r k ^y e y m q k - \ ll]k q^ q k dr (10) 

= p\ + f L(y°p , (y°) 2 q ) + 0(r 2 ) dr . (11) 
Jo 

Here L is a linear function of degree 4 in its arguments with bounded coefficients 
depending on x(t) G M. d . For p the term corresponding to L vanishes except for 
the uniform error 0(t 3 ). Inserting the estimates given in Proposition 13.11 in the 
integrand gives the desired result. For the function Ag;(t) = q^t) — = q + 
tAq l + j Q (t — r)Ag*(r) <ir we set off analogously: 

qi( T ) = G^p J -T^y J ,q e mp a vticn\ a vqi ) (T) = G^p° J -T^ylq e 

gi(r) = ^ Xy > + Gp >-{^ Xy - q )-{^ y - q )-{-^ y - q ) 

Direct calculation gives a long expression the structure of which becomes clear by 
inserting the notation in the proposition. 

q\r) = (Ay)p + {A'y 2 )q + A"q 2 = (Ay)p° + (A'y 2 )q + A"q 2 + 0{t) . 
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A direct consequence is that g*(0) and g*(0) vanish except for an error of order 0{t). 
We find using the estimates of Proposition 13.11 under the integral 

Aq\t) = q^ + tL(y q ,p ) + j\t-r)L(y(T)p(r),y 2 (T)q(T),q 2 (T))dT 

= ql + tL{y»q,y)+eL{yY,{y G ?%,ql) + 0{e) (12) 

where we proceeded as above. In order to achieve an analogous representation for 
the function y(t) the starting equation is given by 

A Vi (t) = Vi {t) - m = Vi + tAy* + \t 2 Ay- + J\t - rfA{ yi )^{r) dr . (13) 
The first order derivatives are explicitly given in (J9j) . For 2y\(t) we have 

{-tf k G km y m y 3 + if (T km y k y m - g jm q m )) Vl - ln^ m y m q l + 9ij {&% - Y^y m q l ) 
= {n"y 2 )% - (n'y) ie q* + 9l] & l Vl = L(y 3 , yq, p) 

with (fi"(?/o) 2 )J = -Ti £ k G km y^ + r{ m Tf y^y^ The corresponding term for jj(0) 
does not vanish but the terms cancel. Here L is a linear function of degree 3 in its 
arguments with coefficients depending on x(t) G W 1 . In the sequel we make use of 
the fact that we may change the names of the indices. This way we get e.g. the 
equality V\ l gj k q k yi = Ti m gj£q e y m . Furthermore 2yf\t) is of the form 

__M i n G n* yjyiyk _ 2 Y^ G mk y . ym _ Y^K k x n y m y k - VtG mk Vmyk 
+ttx n (T km yi y k + g jtf f)y n + Y\ l {Y k ™x n y m y k + 21*™^ - l ]mk x k q m - g ]m q m )Vi 
+T?{T km y m y k + g ]m q m )m ~ ^G m£ x k y £ qi - lijm {Yfx k y^ + G m \y^ + y$)) 
+lijh x k {G^ Vt - T{ m y m q l ) + 9ij (rf x k Pl + & l p t - T^x k y m q e - T{ m {y m q" + y m q")) 

which constitutes a linear function L(y 4 ,y 2 q,q 2 ,yp) of degree 4 in its arguments 
with coefficients depending on x{t) G M. d . For the difference corresponding to y A 
and y A we give a direct method. Expanding the tensor with respect to x and x and 
applying that we are working in Normal coordinates we are left with an error of the 
order 0(t 5 ) for the difference of these terms. For the integrand the estimates given 
in Proposition 13.11 apply. We find 

y(t) = y(t) + tg(x)q + h 2 [g{x)G{x)p° + (Qy)q ] . 

As can be read from the preceding equation the higher order terms corresponding 
to y[t), y\ (t) adequately become part of the error term 5. In order to achieve an 
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analogous representation for the function Ax = x — x we start from the equation 
Az*(t) = 4 + tAx + ±t 2 Axl + f*(t - r) 3 Ai i(4) (r)dr. Here the values of the 
derivatives are determined by inserting into 

x\t) = r%G me y eyj -G%^L(y 2 ,q) 



x 



(3) 



(t) = -(ri n G me + riTf)G nk y k y e y J + riG mi (y e y J + 

= L(y 3 ,yq,p) 



* <( } (t) = - ( + fjjLfsf + f tJn 1 + rJK J G*<y h VKVkvm 

- {YZ n G mi + TZT^)G nk y k (y m + 2y tyj ) + Y^G"' [y, y, + 3y e y 3 + 3y e y J ) 

- G^yf = L(y\y\q 2 ,yp). 

Here L(y 2 , q), L(y 3 , yq,p) and L(y A , y 2 q, q 2 , yp) stand for different linear functions of 
degree 2, 3 and 4, respectively. The explicit expressions of these functions follow by 
inserting the previous results on x, x, and For the difference of the homogenous 
terms x 4 y 4 we proceed as above leaving a uniform error in time in the integrand . 
All together we find 

x\t) = & - -G ij g ijP ° + - h—^G kn G m£ + -^G kn G m n G jn y n y iP ° n . (14) 

This concludes the proof of 13.21 □ 

Using the expression for the solution of the Hamilton equation given in Proposition 
13.21 above we find the approximations for the corresponding derivatives w.r. to the 
initial data go and p° by plain derivation. 

Corollary 3.1. Under the assumptions of \3.1\ the asumptotic behaviour of the partial 
derivatives of the components of the solution to (fjj) is given by 

ap u oqo t op u oq \ox J 

The bounds for the derivatives given in Proposition 13.21 (and hence in its Corollary ) 
can be generalized to higher order derivatives in p° and q which will be used later. 
Let us start with the ^/-component of the solution to the Hamilton equation. Instead 
of the second order Taylor polynomial (11311 we start off with a fifth order polynomial 
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since the second order derivative w.r. to p° needs to be carried out explicitly. This 
corresponds to y^ being a linear function L with argument of degree 6. In the 
proof of Proposition 13.21 the derivatives up to 3 rd order are calculated explicitly. 
Furthermore, up to 4 th order the argument of the n th order derivative of x and the 
{n — l) st order derivative of y are of the same degree. 

For the higher order derivative we proceed in the following formal way, in particular, 
our notation does not distinguish different linear functions L: 

y®(t) = jL{y\y\q\y P ) 

= ^(y 4 ,y 2 q,q 2 ,yp)i + ^(y 3 ,yq,p)y + ^(y 2 ,q)<i + ^r(y)p 

ox ay oq op 

= L(y 5 , y 3 q, yq 2 , y 2 p) + L(y 5 , y 3 q, yq 2 , y 2 p, op) 

Hy 3 q, yq 2 , y 2 p, qp) + L(y 3 q, yq 2 , y 2 p) 
= L(y 5 ,y 3 q,yq 2 ,y 2 p,qp) (15) 

and analogously we find 



y {6) {t) = L(y 6 ,y\y 2 q 2 ,q 3 ,y 3 p,yqp,p 2 ) (16) 

x {6) (t) = L(y 5 ,y 3 q,yq 2 ,y 2 p,qp) (17) 
x {7 \t) = L{y\y i q,y 2 q 2 1 q z ) y*p 1 yqp 1 p 2 ) . (18) 

The algebraic structure behind the procedure is similar to the one studied in |HL01| . 
[NSW85], [BA89J and was developed in the case of a regular diffusion in [KolOOj. 
Since the coefficients are linear combinations of different orders of derivatives of the 
functions g and G we may conclude that the coefficients do not vanish simultane- 
ously. Expanding q and p in a Taylor series w.r. to t we find those functions L which 
explicitly depend on (go) 2 , qoP°, and (p ) 2 . We are interested in the lowest order in 
t for which these arguments appear. 

As was done in the course of the proof to Corollary 13. II for the first order derivatives 
the coefficients of the second order derivatives in q and p° of X and Y are bounded 
given by polynomials in t. This result is summarized in the following corollary. 

Corollary 3.2. Under the assumptions o f\3.1\ the asumptotic behaviour of the second 
order partial derivatives of the components of the solution to (G|) is given by 



> 



&Y = 0{#) 7&?Y = 0{t*) ^Y = 0{t*) 



The result proven in the corollary above can be summarized in the following way: 



d(X,YY 



iK 1 - . . . dC 1 " 
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where we have / = for X and I = 1 for Y as well as Ik = for q and Ik = 1 for 

p° cf. |Koinn| . 

From Corollary 13.21 follows that the matrix JctS is non degenerate hence invertible 
and bounded for each t G (0,t ], < oo. In particular 

d(X Y) 

J(t,x,y,x ,y°) = det - (t,x,y,p,q) p=p o^ x>yiXO>y0) (20) 

U\Pi HI q=qo(t,x,y,x ,y°) 

does not vanish provided the conditions of Proposition 2.1 hold. 



Theorem 3.1. Let the assumptions of Proposition fXIl be satisfied, i) There exist 
positive real numbers c and t (uniformly in xq ) such that for all t < t and \y°\ < 
the mapping (qo,p°) 1 — > (X,Y)(t,x ,y°,q ,p°) defined on the polydisk B c 2 x B E 



c 

t 

3 is 

73" 



a dijfeomorphism onto its image. 

ii) There exists r > such that by reducing c andt if necessary, one can assume that 
the image of the dijfeomorphism described in i) contains a poly disc B r (x) x Bz{y). 

Proof, i) Due to the asymptotic expansion for (X, Y) given in Proposition 13.21 the 
Jacobian of the mapping (qo,p°) (X,Y)(t,Xo,y°,qo,p°) does not vanish for t G 
[0, to]j to > 0, while the assumptions of Proposition 13.21 hold. Hence the mapping 
itself is invertible. Moreover, the mapping constitutes a diffeomorphism since we are 
implicitly given the expansion for the inverse where the t-dependence of the poly 
disk compensates for the different behaviour in time of the parameters. This finishes 
the proof of the statement i). 

ii) We have to prove that for arbitrary (x,y) G B r (x) x Bz(jj) there exists (qo,p°) G 

B c 2 x B c s such that (x,y) = (X,Y)(t,xo,y ,p°,qo). This is equivalent to proving 

7^ 7? 

the existence of a fixed point as for the non-degenerate case treated in |Kol00| . □ 

The projections on the (X, y)-space of the solutions to ([9]) are called characteristics 
of the Hamiltonian H given by (J2J). Suppose that for each initial value (xq, y°) G M 2n , 
and t G [0, to], t > 0, there exists a neighbourhood Q = Q(t,x ,y°) of the origin 
in the (g,p)-space M 2n such that the mapping (qo,P°) (X,Y)(t,x ,y°,q ,p°) is a 
diffeomorphism from f2 onto its image containing a neighbourhood D(xo,y°), then 
the family 

T(x ,y°) = {(X,Y)(t,x ,y ,q ,p°)\(q 0l p ) G Q(x , y )}, < t < t 

of solutions is called a field of characteristics starting at (xo, y°). Due to Corollarv l3.ll 
and the subsequent arguments this means that under the assumptions of Proposition 
I3.1l there exists a field of characteristics for each initial point (xq, y°) i.e. there exists 
a smooth function 

(Qo,P )(t,x,y,xo,y°) : (0,t ] x D(x 0l y°) - 
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such that 

(X, Y) (t, x , y°, (P°, Q ) (t, x, y, x , y )) = (x, y) . (21) 

Recalling that there exists a field of characteristics for any two points (xi,yx) and 
( x 2> 2/2) joint by a characteristic curve of the Hamiltonian H we introduce the func- 
tion 

a(t,x ,y°,q ,p ) = J* ((P, Q), (X, F)) (r) — P(X, Y, P, Q)(r)dr , (22) 

where the components X, Y, P, Q of the solution to the Hamiltonian system corre- 
sponding to H are being evaluated at an intermediate time r, < r < t. Due to 
Proposition 13.11 This function is diffeomorphic equivalent to 

S(t,x,y,x ,y )= f '(Q, P)(r)rf 7 - f H(>y, P,Q){r)dr , (23) 
Jo Jo 

which is called the two point function. 



Proposition 3.3. Under the assumptions of Proposition [Q and Theorem \3.1\ the 

function S(t,x,y,xo,y°) satisfies the Hamilton Jacobi equation 

^ + H (x y — — ) - 
dt ' ' dx ' dy 

as a function of t G (0, to] and x,y G D(xo,y°) with D(xo,y°) as above. Moreover, 
we have 

^( X ,V) = P(t,x,y), ^-(t,x,y) = q(t,x,y). 
° S (x, y), = -p°(t, x, y), ^L(t, x, y) = -q (t, x, y) 



dxo ' dy c 

dqo 



Idea of proof: The partial derivatives are considered first. Since the inverses dx 



and 1 exist due to Corollary 13.11 the functions S and a are identical up to 
the diffeomorphism given in Theorem 13.11 Now the derivatives follow by direct 
calculation. That means that from this point the proof follows the line of arguments 
for the nondegenerate variational principle, see e.g. [iKolOOj. 



Proposition 3.4. Let the assumptions of Proposition HO be satisfied, then the fol- 
lowing representations hold: 



d 2 S 


9P | 


( dX 


dx 2 


dp 1 


[dp 


d 2 S 


9Q . 


( dX 


dxdy 


dp 1 


{ dp 



d 2 S _ dQ ( <9Y N ' 



dy 2 " dq \9q 
d 2 S _ dP ( <9Y N ~ ' 



' dydx dqo V dqo 

where the functions S,X,P,Q depend on the same tuple (t,x,y,qo,p°) of variables. 
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This is achieved by iteration from Proposition 13.31 and Corollary 13.11 which implies 
Theorem 13.11 

Remark 3.1. Due to the symmetry of the results in Corollary 1 3. 1\ and Proposition 
\3.3\ an analogue to Proposition \3.4\ in the variables Xq, y° holds. For the lowest order 
approximation the following explicit expressions exist for each component: 

P{x,y}-\=£(1 + S) 1(1 + 6)) ^{x ,y°}"U( 1 + 5 ) 1(1 + 6) J - [ ) 

For the Hamiltonian H(£,r)) let us define the functional 

J t (0 = J max { ( V , - #(£, V ) } (r)dr (25) 

The maximum is taken over all piecewise smooth curves £ with given endpoints 
£(0) = £ and £(t) = x. 

Proposition 3.5. For the Hamiltonian given in p|) the characteristic curve 

(X, Y)(t, xq, y°, qo,p°) provides a minimum for the functional ( f%5)J over all curves in 

D(xo,y°) = B c 2 (q ) x B c a(p°). Furthermore, S(t,x,y,xo,y°), of {23}) . constitutes 

the unique minimal value of the functional §2i 



Proof. For the Hamiltonian given in j2]) the Weierstrafi condition is satisfied, i.e. for 
all t > and all £, r) e R d : W(x(t),y(t),£,rj,p(t),q(t)) > for arbitrary solutions 
(x(t),y(t),q(t),p(t)) of ((9)) where the Weierstraft function is given by 

W(x,y,Z,r],p,q) = AH-L(H) 

with H as defined in (J2J), AH := H(x,y,qx,pi) — H(x,y,q2,P2) and L(H) = 

( A <?> f ) + f ) : = (ft - f ) + {Pi ~ ft. f ) • T ^n we find 

A# = ~ (<7(a:)gi, gi) - (#(x)g 2 , q%) - (G(x)y, Ap) - f-^-(G(x)y, y), Aq 

L(H) = (q x - q 2 , g(x)q 2 ) + ~ (J-(G(x)y, y), Aq^j - (G(x)y, Ap) . 
Combining the two equations reveals 

W = ^(9(x)qx,qi)-(qi,g(x)q 2 )+-(g(x)q 2 ,q 2 ) = ^ (qi-qz, g(x)(qx-q 2 ))+0(l) > 0, 

the positivity holding since the matrix g is positive definite. As a consequence the 
Hamilton function treated in this paper satisfies the Weierstrass condition. Inserting 
definitions and Theorem 13.11 we have If > S c.f. (1231 . Hence the minimum for the 
functional will be unique for arbitrarily chosen characteristics in B c 2 (q ) x B c z (p°), 

which finishes the proof of Proposition 13.51 □ 



12 



Combining the result of Theorem 13.11 with the uniqueness result obtained in Proposi- 
tion [33] gives us the equivalence of solutions to the Hamilton equation and Hamilton 
Jacobi equation. 

We conclude this section by showing that the two-point function S is convex on 
D(xq, y°). For arbitrary square matrices a, b, c, d the general following formula holds 



Using the asymptotics given in Remark IBTTl as a starting point, applying the above 
result and integrating twice w.r. to x and y we see that the two point function can 
be uniformly approximated by a quadratic function in x and y. In particular the 
two point function S can be approximated as follows 



This leads to a Gaussian approximation of the solution to jTJ). 

4 Uniform Estimates for the Heat Kernel 

In this section we shall construct an asymptotic expansion for the fundamental 
solution (Greens function) of the degenerate diffusion (pQ). To this end the well 
known WKB-type method is modified in such a way that the variables representing 
derivatives i time of order zero to three are being transformed in order to lift them 
to the same time scale. Hereby closed algebraic equations for the coefficients of the 
expansion are achieved, see also the construction in [KolOOj for the origin of our 
construction. 

We shall discuss two types of asymptotics for the solution to (JTJ) and (J2J), namely, 
small time asymptotics when t — > and h is fixed (h = 1) and small diffusion 
asymptotics when t is fixed and h — > 0. 

We look for the asymptotic fundamental solution to (J3| using a WKB Ansatz 

u°(t, x, y, x , y°, h) = C(h)(j)(t, x, y, x , y°) exp[-S(t, x, y, x , y°)/h] , (27) 

with real functions and S called the amplitude and the phase, respectively. More- 
over, C is a function only depending on h alone while (f> and S are independent of h 
having the arguments x,y,x ,y° G M d . This particular form of the asymptotic fun- 
damental solution applies due to the positivity of the resolving diffusion operator. 
It is specific to linear second order pseudo differential equations describing Markov 




det (a — bd c) . 



S(x,y,x ,y ) 



6 



(1 + 6){x- x ) 2 + -(1 + 6) ((x - x ), y + y°) 
+ hl + S)[( y , y ) + ( y , y ) + ( y °, y )]. 



(26) 



processes. 
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This asymptotic ansatz is also called exponential. Hence one should require that 
the solution u of (OQ) differs from u° by a multiplicative error: 

u(t,x,y,x ,y°, h) = C(h)<fi(t,x,y,xo,y )exp[-S(t,x,y,x ,y )/h](l + 0(h)). (28) 

For (1251) only the term with minimal entropy of the sum survives at each point. 
Therefore, for ((5|) the common superposition principle transforms to the idempotent 
superposition principle (Si, Sz) min(Si, S2) at the level of actions. For a detailed 
discussion of this idempotent superposition principle and its applications see [KM94J, 
|KM97j . Inserting u° ((271) into (OQ) reveals 



, (d<j> I ± dS\ h 2 , Jd 2 $ (j)d 2 S\ L N dcj) 1 ± dS\ fnn . 
-h ( d , . d</> 1 A dS\ 1 / . . dS\ , / . . dS d<p 

The terms in the above equation may be classified according to the degree of h which 
ranges from to 2. Comparing coefficients we obtain the following two equations: 
The first Hamilton- Jacobi equation (h to the order 0): 

+ o 9(x)^, -7r)-\ G ( x )v, IT + o -^-(G(x)y, y),— )= (30) 



<9i 2 \ dy 1 dy J \ ' dy J 2 \dx ' 1 dy J 

which in terms of the Hamiltonian function ([2]) gives 

d a TT f dS dS\ n 
at \ ax ay ) 

This way the system given by the PDE P| is linked with the Hamilton- Jacobi equa- 
tion studied in Section 2. The second equation is the so called transport equation 
(the coefficient of h 1 ): 

(31) 

Going the opposite way by choosing S and which satisfy (130|) and (l3il and inserting 
then into the r.h.s of (OQ) we define a function u° which satisfies equation (pQ) up to an 
error of order h 2 . More precisely we obtain the following inhomogeneous equation 

du° f d d\ h 2 

h—-H[x, y, -h—, -h—J u °=--C(h)F (t, x, y, x Q , y°) (32) 

F (t, x, y, xq, y°) = tr{q(x)^^\ exp(— s< ^ t ^ y ' > ) i. e . the right hand side defines a func- 
tion F (t, x, y, xq, y°). If we would concentrate on a Cauchy problem instead with a 
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smooth initial function then we would have to solve an appropriate Cauchy problem 
for the Hamilton-Jacobi equation. 

We now proceed to derive the precise asymptotic behaviour of the solution to the 
parabolic equation (0Q) for small times (respectively in a small neighbourhood of the 
starting point) we shall provide asymptotic expansions in the time parameter t for 
the two point function S and the amplitude given in (1271) . Studying the inho- 
mogeneous system (l32l) has the advantage that the functions S and </> satisfy the 
Hamilton-Jacobi- fl30|) respectively the transport equation (IBTl) . 
In Section 3 of this paper we have seen that the two point function S in [23] is 
smooth and satisfies the Hamilton-Jacobi equation almost everywhere. Moreover, 
it was shown that the solution of the Hamilton-Jacobi equation and the Hamilton 
equations are equivalent thus generalizing the results given in [KolOO] for non de- 
generate systems to a special case of a regular Hamiltonian of rank one. Due to the 
results in Section 3 we only use the expansions of the solution to the Hamiltonian 
system in the sequel. 

In a fist step we neglect the inhomogeneity in (132]) . Starting with the phase respec- 
tively the two point function we are looking for a representation of the form 

oo 

S{t,x,y) = Y,Si(x,y)t i . (33) 

i=— k 

Traditionally this is achieved by introducing the asymptotic expansion for the cor- 
responding Hamiltonian system j9]) into the Hamilton-Jacobi equation at the point 
(t, x, y) and subsequently solving a closed system of algebraic equations in x and 
y G M. d . This fails in the case of our degenerate system. In this paper we proceed 
by introducing a change of variables which first shifts the origin into the solution 
(x, y,p, q) of the Hamilton equation ([9]) starting at (x , yo, 0, 0) and then adjusts the 
time scales of the variables on the manifold and the tangent space. We define a 
function ^onl + xR d x R d by 

Y,(t, Z, V) = S(t, t(x + x(t)),y + y, x , y°) . (34) 

We remark that the function a(t,x,y) = S(t,x + x(t),y + y,x ,y°) coincides with 
the one given in ([23|). There holds ^ = and f| = ^ + . In this way a 
closed system of algebraic equations for the coefficient functions Si of an expansion 
of the phase function in powers of x and y will be derived. 
The equation (1301) becomes 



° = ~dt i d£ - [ 9{x)q > dy- ) (35) 



1 

+ 2 



-(y + y),(y + y)) - —^-ViV 



dx ) V dx 



as 1 / , , ,dY> dY\ 
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Proposition 4.1. For the transformed two-point function E as defined in (3$ ) there 
exists a unique asymptotic expansion in powers of t > of the form 



t 



S + t£i + t 2 E 2 + 



such that Ei anc? E vanish at the origin, Si «s strictly convex in a neighbourhood 
of the origin and all Sy, j > 2, are regular power series in (£,?/). Moreover, the re- 
mainder of arbitrary order can be estimated by means of induction using Proposition 
EH as initialization. 



Proof As in Section 3 let us replace the coefficient functions g, G in equation [35] by 
their n th order Taylor series and let us also insert the n' th order approximation of 
(x, y,p, q). Then the asymptotic result given in Corollarv l3.2l gives that the combined 
error due to this replacement is of the order 0(t n+1 ). 

This way (1351) becomes in coordinate form (normal coordinates) for n = 1 up to an 
error of order 0(t 2 ): 

^ (g + y) % + fggfe - yl) - yg) ( Vj + yj) - y° k y Q e y]\ + Q(t 3 ) gS 
dt t d& 

- ^tgfM' - VIM + »?)(% + »°) - + °< i2 )]^ < 36 > 

+ ;d + >-*(e.-rf) + o(<f = o. 

We now rearrange the l.h. side of this equation according to the order in t. For 
the lowest orders in t we arrive at the subsequent set of equations. In particular 
for the lowest order in t, i.e. t~ 2 , equation (l36l) reduces to the following first order 
homogeneous partial differential equation with linear coefficient functions: 

~ ~ {y + + 2 {9 °^y-' ~W ) = Q - (37) 

This linear first order partial differential equation is of the following general form 

Am + (Ax, Vu) = p(x) (38) 

where A G IR + , x G IR m , p(x) is a homogeneous polynomial, and A is an mxm matrix 
with strictly positive eigenvalues a 1 < a 2 < . . . < a m . The solution constitutes a 
polynomial of a particular form given in the following lemma. The proof goes by 
inserting and direct calculations. 
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Lemma 4.1. i) Let p be a homogeneous polynomial of degree q. Then the solution 
of exists, moreover it constitutes a polynomial of degree q with coefficients of 
the form 



d " u - 1 (c-x )... (c-i ) dq p c C - 

dx h . . . dx lq a n + ... + a jq ™> {U '* } dx tl . . . dx iq " ^ 



Let p = Y^=oPg be a sum of homogeneous polynomials of degree q (in case m = oo 
the sum has to be absolutely convergent in a ball). Then the analytic solution of 
(ESP exists and is given by the sum Y^=i u q °f solutions u q corresponding to the 
inhomogeneity p q . Ln case m = oo the domains of convergence of p and the solution 
u to [38\) coincide. 



For a proof see [KolOOj. We now proceed to calculate the coefficient functions Sj, 
-1 < i. 

Lemma 4.2. Under the additional assumption E_i(0, 0) = and that E_x is strictly 
convex at the origin the formal power series for S_i in powers of £ and y reduces to 
the quadratic form 6(£, £) + 6(£, y) + 2(y, y). 

Proof. Applying the assumption that £_i has an asymptotic expansion 

fc>0 

where (k is a polynomial of order k in ^ and y and using the assumption £_i(0, 0) = 
it immediately follows that Co — 0. In local coordinates g — 1. The first order term 
is of the form £i = + Ly where K, L 6 R. By inserting into (1371) and using that 
Co = we get the equation 

-K^-Ly-(y + OK + h 2 = 0. 

At the origin we find L 2 = 0. Differentiating with respect to y reveals —K — L = 
and hence £i vanishes altogether, i.e = 0. Assuming for the quadratic part (2 the 
representation 

-l(A£,0-(Bt,y) + ~(Cy,y) (39) 

we get the equations 

A = -B f B 2B + A = CB -C + B = -C 2 . 
2 2 2 



Solving this algebraic system of equations we find, since (13911 should be positive 
definite, C — 4 id, B = 6 id, and A = 12 id, i.e. (1391) is of the form 



6(£,0 + 6(£,y) + 2(y,j/) . (40) 
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Finally all k > 3, vanish, in particular for = 3 the terms of order 3 in £, y make 
up the equation 

o . _ Cs _ ( (y + , |) + (_* _|) . C3+ ( (9 + G , | _ « + 4s) , | 

(41) 

where f[40l) has been inserted and the right hand side has been multiplied by —1. 
Due to Lemma |4~T1 this implies that ( 3 = 0. Inserting £fc, fc > 3, into f|37l) reproduces 
(JlTjl . hence Cfc vanish f for k > 3, which finishes the proof of the lemma. □ 

For the term S corresponding to the order t° in the expansion of S we have the 
equation 

-to+f)^ + («e+^ = o, 



the solution of which vanishes due to Lemma l4~Tl since the equation is homogeneous. 
For the coefficients of t in (1361) we obtain the equation 

^i-(y + 0^ + m + ^y)^- = m,y) , (42) 

where Fi(£,y) is the sum of the homogeneous polynomials Fn,Fi 2 and F13 of order 
2,3 and 4, respectively, in £ and y. In particular we have 

F n = 4?[(12&& - 4y iyk )y?y¥ - (18£ <% - + 7 W + 9^)^%° + (3£*6 + 2&y,)l/M] 

F12 = $[(-6£i£k£i + %k£ilVi + 4£kyiyi)i/j 

+fi6€i€kVj + iHkViVj - tyiVjVk + 1 )//J r 
^i3 = gS^hViVjyi - 4£k£iyiyj - 18&&6%- - 9&&&6] . 

Due to Lemma I4TT1 the form of the solution £1 is determined to be a polynomial of 
degree 4 in £, y and y°. The expansion of £ in powers of t (having expanded in x 
and y before) is explicitly given by 



+ 1 \-jE9i,jyo,kyo,iyiVj - 20 91 j vojeVo^iVj + 5 9ijyo,kyo,i£itj 

- s 9ijyo,jyo,iyiyk + 1 gtjyojyo&yk + f stjyojyo^k + j\ gijyo^ojyiyk 
+ 20 9i,jyo,iyoAkyi + s&j-iAvvoj&sj + to 9ijVo&tVkyj - To 9i,jyo,i&€jyk 

3 k.l r r r 1 k,l (- 6 k,l s- f- feJ 3 Mr* 

- 2 9i,jyoA£j£k - s 9ijyojQiViyk - 5 g^yoj^kyi -^g^yoj^A - ^g^iiyiykyj 
+ 1 9i$U s vm + f <#&&&yj + + 0(t 2 ) . 

The calculations for higher orders go by analogy although they become substantially 
wilder. The terms £j, j > 1, turn out to be homogeneous polynomials in £, y, y° of 
degree j + 3 which can be estimated analogously to the term Si. 
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In normal coordinates in particular the term Si is given by 2.5#^'y°20°y° due to 
the symmetry of which is why we omit it when transforming back. The inverse 

transformation S(t, x,y,y°) = x ~ x . ® , y — y(t)) has been found explicitly for an 
approximation of order 1 using MAPLE, a programme for symbolic and numerical 
algebra, as being 



Si(t, x, y, x , y°) = t' 1 (-6 ^ - 6 + 6 % -f + 6 
- 12^ + 6^ + 6^ + 2^ + 2^ + 2^ + 0(0 



This is consistent with the estimate found in (l24j) . On the diagonal this expression 
is of the form Si = Gyft' 1 which obviously is not vanishing and hence reveals an 
exponential decay for the kernel on the diagonal. 
For Xq = this reads: 

S(t,x,y,y°) = ^[l + 0(\x\ 2 ) + 0(\y-y°\ 2 )]\x\ 2 
+ ^(x,y + y° + 0(\x\ 2 )) 

+ - t l\y\ 2 + (y°, y) + |y°| a + 0{\x\ 2 ) +0(\y- y°\ 2 )} ■ 

Having completed the first step in giving an approximate solution of exponential 
WKB type (|27|) . namely, having derived an asymptotic expansion for the solution 
5* to ([301) we now turn to the transport equation (I3T1) . The parabolic equation fl3Tl) 
has the amplitude (j) as the only unknown function left, for which we intend to give 
an asymptotic expansion in the sequel. As before we introduce a change of variables 
which adjusts the different time scales for the variables of the tangent bundle, in 
particular: 

£, y) = tf + x(t),y + y(t)) (43) 

where the constant a is determined by consistency requirements. In terms of the 
transformed function the transport equation takes the form 

= i - t* + 2 tr +x) w) +(m+x) ww- <s(xK w (44) 

Bib 1 r)Q Bib 
-r\t + G{ti + x)(y + y)- G(x)y, ^) + + x)(y + y), (y + 

Proposition 4.2. There exists a unique a > smc/j i/iai t/zere exists an asymptotic 
expansion for the solution of ftBl 



1p = l/l + tip! + + • • • 
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where each tpk, k > 0, constitutes a polynomial in (£,?/) with degree strictly smaller 
than k and ip® is some constant. Moreover, the solution is unique up to a constant 
multiplier. 

Proof. Assuming the solution of (1441) to be of the form of an asymptotic expansion 
(regular power series) the equation (jUJ) decomposes into separate equations for each 
order in t, e.g. for the lowest order (t -1 ) we find 

1 <9 2 £ 
-aipo + -ipo tr(g(x ) Q ~ {x ,y )) = , 



where we also have inserted the asymptotic expansion given by Proposition ^. 11 This 
determines the parameter a in fl43l) to be 



1 <9 2 £ 

" = j tr (9( x o) d 1 {xo,Vo)) ■ (45) 

Recalling that we are working in normal coordinates and reminding the representa- 
tion given for £_i in Lemma B~2l we find a = In. Comparing the coefficients of the 
term t° reveals the following equation 



the solution of which vanishes due to Lemma 14.11 In particular, this equation is 
the homogeneous counterpart to ([321) . For the coefficients of the term t 1 we get the 
following nonhomogeneous counterpart of the previous equation: 

* - (i + v, f ) + (ee + %, f ) + tr (I^fi + - rf) (6 - rf)) = o. 

Decomposing the trace into its homogeneous polynomials of degree 0,1,2, respec- 
tively, we may determine the solution (a polynomial of degree 2 in (£,y,y )) as 
indicated by Lemma I4~T1 

^ = 1 + e (- (-! &k - 1 9 J k %k - 1 & + 1 rfjyoj - 1 .</;;>),, 

3 M 3 iA 3 i 5 i fc.Z \ , 5 / 7 i.fc 1 i.i 3 £,?' 

- 56 9k,iVj ~ 56 - 56 9kjVj ~ 9i4 Vo,l )Vk + Q [jo g-jVoj - 5 9kjy0,j - 56 9k,iVj 

k,l 



3 1,i 3 i,i k,l \ r 1 fc.i 5 k,l t- . l «,fe 

56 _ 56 ffe,j%- _ S/0,1 J £fc - 6 9i,i VofiVl ~ g !/0,fc6 + 210 9ijVkVj 

[40 & j W + 35 9i,j ZjZk + go 9i'i VWk + 20 ZkVl + 5 4fc& 

k9i'iVo,kyo,i + T5 9ijyo,jyo,k - ^ ^yo^yoj) + 0(t 3 ) . 



(46) 
□ 
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The inverse transformation <p(t, x, y) = ip(t, 



-x{t) 



y — y(t)) has been done explicitly 



for an approximation of order 1 using MAPLE. Applying this back transformation 
to if) (which is given by an asymptotic expansion in t) will give back an asymptotic 
expansion in t. In particular, for x = we have 

<f>{t, x, y, y Q ) = (1 + 0(\x\ + \y - y°\ + t))t' a . 



Inserting the expansions for S(t, x, y, Xo, yo) and into the WKB Ansatz (1271) reveals 
directly that the leading term is Gaussian. 

We now proceed to show that the error term can be expressed by a factor. For 
an arbitrary m-dimensional system of ordinary differential equations -^X(t,a) = 
f(X, a) depending on an m-dimensional parameter a a result of Liouville states 
that the square matrix J = [^]y is non degenerate on some time interval and 
satisfies j- t J = J £f[f^]- For a proof see e.g. |FM81| . It is possible to determine 
the amplitude explicitly by applying this result to the system 



d dH / 



and 



where H is given in (J2J) and by using the equivalence of solutions for the Hamilton 
and the Hamilton- Jacobi equations given by Propositions 13.31 and I3.2L We find 

P = ^;(t,x } y,x 0} y ) and q = ^-(t, x, y, x , y°) 

to be dependent variables of x, y. Hence along the characteristics (X, Y)(t, x , y°, q°,Po) 
of the Hamilton equation corresponding to the Hamiltonian ((2D the determinant 
J = det g^'p^ satisfies the equation 



J = Jtr 



d 2 H d 2 H 



dpdx 
d 2 H 



dpdy 
d 2 H 



dqdx dqdy 



+ 



d 2 H d 2 H 

dp 2 dpdq 

d H d 2 H 

dqdp dq 2 



d 2 S 
dx 2 
d 2 S 



d 2 S 



dxdy 


dydx dy 2 



where the second term in this sum corresponds to the dependent variables p, q. 
By inserting the definition of the Hamiltonian (J2J) most partial derivatives of H 
vanish leaving us with 



J — Jtr I — 



dG 



dG d 2 Hd 2 S 



dx ^ ^ dx ^ ^ dq 2 dy 2 



Jtr I g(x 



d 2 S 
dy 2 



(47) 



which implies 



d_ 

It 



J 



1 T _i 
■-J-* tr 



d 2 S 
dy 2 



(48) 
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On the other hand inserting the Hamilton equations into the transport equation 
f|3Tl) reveals 



d 

di' 



+ 



1 



tr 



d 2 H 

dq 2 



x,y, 



dS dS\ d 2 S 



or more explicitly 



dx dy J d 2 (x, y) 
d 2 S 











(49) 



Comparing (l4"8l) and (I4U1) we find that and are identical up to a constant 
factor a. 

We now deal with the inhomogeneity in the PDE (1321) . Aiming at a multiplicative 
error for an exponential approximation of WKB type respectively a multiplicative 
error for u° we need uniform estimates for the inhomogeneity in an appropriate 
neighbourhood of the starting point (x , y°). This will be achieved by exchanging 4> 
with J~\ and using the estimates given in Corollary 13.11 and ffl~2)) . 



Theorem 4.1. Let us assume that the assumptions of Theorem lff.il are satisfied 
in the whole cotangent bundle. Let an approximating solution of the form 
to (Qy satisfy the Hamilton Jacobi- (Sty) and transport equation ( EZP - Then the 
inhomogenity F fffgj) is a multiplicative error of the type 



h— - H [x,y, -h—, -h— 
at \ ox ay 



2„,0 



M u = 0(h 2 )t 2 u 



Proof. Recalling the definition of F and the fact that and J 2 are identical up 
to a constant factor c.f. (T48l) and (l49l it suffices to examine 



-J z(t,x,y,x ,y°) and 



-J i(t,x,y,x ,y°) 



[dzt)"" > (d Vl 

for v = 1, 2 and 1 < i < n. Let p = p {t,x,y,x ,y°) and q = q (t,x,y,x ,y°) 
satisfy the Hamilton Jacobi equation (1301 . Then we have along the characteristics 
(X, Y)(t, Xq, y°, q°,Po) of the Hamiltonian given in (J2J) for 7 = (70,71) with 70 := x, 
71 := y both in R d and ( = (q°,Po) G 

d 



p2<2 



8% 



■J-5 



1 



1 



d 2 (x,r) 

d(d( 



1 in 



■J 2 tr 



y) 



-1 



djX, Y) 

dC 



(50) 



since 



d(X,Y) 



H 2 " 2 57, 



9 d'c ' Analogously, by iteration 
1 9 



0% 



3 T 5 . d tN2 T 3 d 2 

-J * J f + J 3 

2 ^ #7; 



mfc 



■J 



— J-2/sT 
2 
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where 



K -I Hi TV + r-l ° 2J ^ j-i_d_ T d 2 C l 




+ 



+ 



3 ( d(X,Y) 
2 V dC 



Here the last factor may be rewritten as follows = — 0§f) d( l d( n (ic) ' ^ n 
order to simplify notations we replace I m by m, etc. in the next formula where 



I m was defined in ([19]) . Applying Corrollary 13.11 we find the following asymptotic 
behaviour for K: 

0{t 2 ) + o(r M )o(t M " w -^)0(r (w " i) )o(i _(!H, - 1 ' ) ) + 0(r {3 - k ~ m) ) 
= o(t 2 t) + o(t 2 t) + o{t 2 i) . 

Recalling that % = 1^ = 1 for yi, 1 < i < d, and inserting into (1501) reveals for the 
leading terms 

= j-\ {t l ) (51) 

and by iteration 

^J~\ = J-lO(f) (52) 

which finishes the proof of Theorem 14. 1L □ 

Next we relax the assumption that Theorem 14.11 hold in the whole cotangent bundle 
by introducing a cutoff function in the definition of the function uq. Let Xd be a 
smooth function with values in [0, 1] such that \d vanishes outside of the domain 
D of the cotangent bundle, is 1 for all points (x, y) in the domain D except in a 
neighbourhood of the boundary dD of D, moreover, all partial derivatives of xd 
until order one in x and until order two in y are bounded for t G (0, to), to > 0. 
Then the function 

u D = C{K)xd{x - x ,y -y )(f>{t,x,y,xo,yo)ew[-S{t,x,y,xo,yo)/h] (53) 
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is globally well-defined for t < t and satisfies 



h 



du D 



h 2 

— tr 
2 



,2„,£> 



9{x) 



d u 



h I G(x)y. 



du 



D 



dyidyj 

= -h 2 C(h)F(t,x,y,x ,y°) 
where the inhomogeneity is now given by 



dx 



h 

+ 2 



d 

—(G(x)y,y) 



du D 

dy 



(54) 



tr G(x) 



d 2 



+ tr g(x 



dyidyj 
9xd 



Xd exp 



S 
h 



+ 



tr g(x) 



dyj dyi 



+ tr g(x 



+ h 



-2 



G(x)y. 



D 



dx 



| 9 2 xd 
dyjdyi 



1 / d 



exp 

9xd 
dy 



9xd dS 
dyj dy { 
S_ 
h 
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In the sequel we are going to study the order of magnitude in t and h of the inho- 
mogeneity F. First we collect some general assumptions and facts. The appropriate 
neighbourhood is the polydisk D = B r (x ) x Bz (y ) in particular its closure D. The 

functions xd-, ~§^Xd, ~§^Xd and -§^Xd are bounded on D, moreover, g and G = g^ 1 
are bounded since the manifold is compact. In normal coordinates (actually x = 0) 



g(x) = id + GijXiXj + 



. . . and G(x) = id — Gij(x — x )i(x — x )j 
,„.:,. ■ x ) =O(-y t) = cAr and 0(£G) = O (- 

c A r. Using the expansion given for S in (1261) we find 



that 0(£g) 



0(x 



+ . . . which implies 

x )) = O(y t) 



[x 



dS_ 

dyi 



5)(x — xq) 



t 2{ 

6 2 

-(1 - 5)c A r + -(1 - 5){y + tq + y + 0(t))i. 



this implies 0(g) = 0(t- 2 ) + 0(r 2 ; 
similarly 0(g) = O(g ) = 0(r 2 ). 



o(r 



Using Proposition 13.31 we get 
Applying (1501) to J^- and collecting terms 



Theorem 14.11 is modified by the cutoff xd showing the following order of magnitude 
in t and h for the inhomogeneity F: 



0(i 2 )u D + 



< 0(t 



u D + 



+ 0(t 1 )0(r 2 ] 

(0(r 2 ) + 1 + 0(t)) X AD s O [ r << , 



(o(r 2 ) + o(t) + i) + /*- 2 (0(t 
h- 1 





S~ 


^jX ad J exp 





where Xad £ stands for -§^Xd as well as -^j-Xd- The support AD e of these functions 
is a neighbourhood of the boundary 3D of width e. In the above estimate we used 
that for given 5 > 0, h > one can choose t such that for t <t and h < h 



12 n (2nt)- 2n < exp 



_5_ 

hi 3 



ht 3 



< exp 



_5_ 
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The constant Q = Q + 5 summarizes the constants Q = mxn{t 3 S(z, £)| \x — x \ = 
r — e and \y — y \ = and 5. These calculations are summarized in the following 
proposition. 

Proposition 4.3. Under the assumptions of Theorem \3.1\ the inhomogeneity F 
has the asymptotic behaviour 0(t 2 )u D + 0(exp[— j^3])xad e for t, h > and some 
positive constant Q. 

After having studied the solution (1531) of WKB shape we passed on to solving the 
inhomogeneous equation (1541) . We now reconstruct the solution to the original homo- 
geneous equation using the well known DuHamel principle. Recalling the classical 
results of Egorov, Hormander, Oleinik and others we know that equation ([1]) is 
hypoelliptic since the generator £ is of the form 

C = S"=i^ 2 + Xo 

where X t = A tj -^-, 1 < i,j < n with A* A = G, detA ^ 0, and X = (G(x)y, J^) - 

\ ^(G(x)y,y),^y moreover, the Lie-brackets [X u (G(x)y, £)] = (A*)"/^, 1 < 

i < n, together with the vector fields Xj, 1 < i < n, form a basis of the tangent 
space. In particular it is shown cf. Hormander |Hor68| . Oleinik |Qle71| that for any 
initial distribution in D' and smooth inhomogeneity there exists a unique smooth 
solution to (Jll). From this it immediately follows that the operators (e* £ ) f>0 possess 
a smoothing kernel. For the sake of a shorter notation we are using z = (x, y) and 
zq = (x ,y°). 

Proposition 4.4. For any domain D in IR and C 2 function f with support in D 
there exists a unique smooth solution v t in D satisfying the inhomogeneous equation 

9 r t 

In fact, v is given by the integral formula 

v(t,z,z ) =u(t,z,z ) - / u{t-T,z,£)f{T,£,z )d£dT 

JO JM. 2n 

where u satisfies the corresponding homogeneous equation ^J\). 

Proof. Uniqueness follows from the uniqueness of the corresponding homogeneous 
equation ([J]). As for existence, applying the operator — C to the r.h.s. above 
yields 

(I - c ) v = (I - c ) M(t ' z ' zo) - (I - c ) I L u{t ■ r ' *' 0/(r ' e ' zo) didr 
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f f u(t-T,z,£)f(T,£,z )d£dT + £ [ [ u(t-T,Z,g)f(T,Z,Zo)d£dT 
OK Jo JR 2n J0 </R 2 ™ 

u(0,z,€)f(T,€,Zo)d£- [ [ ^-u(t-T,z,£,)f(T,£ } z )d£dT 

Jo Jm. 2 ™ at 



+£ u(t-T,z,£)f(t,£,z )d£,dT 
Jo Jm 2n 

fit, Zq, Zq) . 



We have used that u solves the homogeneous equation to eliminate terms. Further- 
more when differentiating the parameter integral w.r. to t we have skipped the reg- 
ularization needed for the argument giving J R2n S z (£)f(t, £, zq) d£ and have exploited 
the cutoff given by the function / for obtaining the integrability of (J^it)/. □ 

For the homogeneous equation (pQ) and its nonhomogeneous counterpart (1531) the 
solutions are therefore related in the following way 



u D = (1 - hl t )u , 

with u D given by (1531) and where the operator I t is defined by Itu(t, z, z ) := 
Jo J R 2nU(t — r, z, £)F(t, £, Zq) d^dr with F given in (154"1) . In case the inverse op- 
erator exists, which is equivalent to the convergence of its associated formal series 
I + ^2 h m I™, we retrieve the solution to (DQ). 



Theorem 4.2. For small t the fundamental solution u to (CP has the form 

u = u D (l + 0(ht)) + 0(expi-^L]) 
with u D given by p5Jj and Q as in Proposition ^. 3 . 



In order to meet the assumptions of Proposition 13. 11 Theorem 13. 11 and Proposition 
14.31 we choose D = B r (x ) x Bi_(y ) for r sufficiently small. Due to Proposition 
14.31 the action functional S satisfies min^{S'(t — r, z, £) + S(r, (,,z )} = S(t,z,z ). 
Introducing a zero term, using that S is a convex function on D with a unique 
minimum, and replacing S(t — t,z,^) + S(t, £, z ) — mm^{S(t — r, z, £) + S(t, £, z )} 
by its Taylor polynomial of order three at £o minimizing S(t — t,z,£) + S(t,£,zo) 
we find: 

ku D {z)<C [ [ u D (t-T,z,o(r 2 u D (r,^z ) + e-^XAD^-^)) d£dr 

Jo JR 2d V 7 

< ^ S(t ' Z ' Z0) f u {h) ™„ I e~^-^M-^U-^ XD ^ _ Zo)xd(z _ d ^ dT 



JO \ l 1 I ' JMi^ a 

+ C{h) If (t- rr 2d e~r^r, z ,o XD{z _ £) e -^ XADe (£ - z ) d£ dr 
Jo Jm. 2d 



'0 JR 2d 
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< CC(h)r 2d e-^ sit > z ' zo k X D*i {z - zq) 



+C{h) 



-M 



(t _ T y^-i^-0^) XD{z _ t) XD ^ z - ZQ ) d ^ dT 



< C (C(h) 2 e-^ z °h- 2d+1 + te"^) X d*< z - z ) 



(55) 



where C{h) 



12'' 



(iirht) 2 



and Xd* 2 ( z — z o) is the characteristic function of the poly 



disc -D* 2 = {(x, y) G M. 2d \ \x — Xq\ < 2r, \y — yo\ < y}. Moreover, above we have 

' and = * - ^ - d ' having 



introduced the matrices: 



ar',2/') G 



non vanishing 2 x 2-matrices on the diagonal as given in ([2l]) . where £ 
IR M . The determinants of both matrices are given by tj with <r G {£ — r, r}. Their 
sum equals t 2d r~ 2d (t — r)~ 2d . All constants are included in the universal constant 
C. Finally Q is given in Proposition 14.31 We have used Xd(£, — zo)xd{z — £) < 
Xd+ 2 ( z — zo) and that the function exp(— -J^) is isotone having the image [0, 1). The 
m th order term I™ follows by iteration. For the sake of a clear representation we 
introduce the following notations: 



with z 



f{t,z,z 

9m(t, Z, Z 

(x,y) G R 2d , 



t 2d XD(z - z )C(h)exp(-S(t, z, z )/h), 
Xd+™{z - z ) exp - — 

(xo,yo) G M, 2d and Xd*™(z — z ) is the characteristic 



function of the polydisc D* m = {(x,y) G ~R 2d \ \x — xq\ < mr, \y — yo\ < ^f 1 }. In 
terms of the new notation inequality (1551) reads / * f(t,z,z ) < CC(h)f(t,z,z ) 
together with / * Cgi(t, z, Zq) < Ctg 2 (t, z, Zq) for some constant C . As to the higher 
order convolutions we apply: 



\f * 9m{t,Z,Z )\ 



XD+™{z-£)e 

JR 2d T 



Xd{£ ~ z ) d£dr 



< Xd+™+Az~Zo) / e Ht-r)* e - 

Jo 

< g m+1 (t,z,z )e~ M t. 



M f C ( h ) (3 f . t (Z-zoU-zo) 



-2d 



Xd{£ ~ z ) d£dr 



where we have used that S is locally convex satisfying S(t,z,z ) > M + (Hy >t (£ 
Zq),£ — zo) on D. Moreover, there holds 

-t 



91 *9m(t,Z,Z ) 



JR 2d 



Xn*m(z-^)e m*-^ e ^Xd{£~ z )d£dT 



< XD* m + 1 ( z — z )e hi 5 



- n 
e h^ 5 



Xd{£ ~ z )d£dr 



< 7r n r(- + l) r d XD * m+ i(z - Zo)e~^ 



e'l^dr < tb(r)g m+1 (t, z, z ). 
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n n n-s 



where we have estimated t n e f^ 5 < t n e < e "w 3 " < 1 and b{r) summarizes 
the constants. Hence by induction for the m th order term we find 

/i m - 1 /® m " 1 u D (z) < CC(h)(Chb(r)t) m ~\f + 2g 2 + . . . 2 m ~ 1 m ) 
< CC(h)(Chb(r)t) m ~\f + g n 



This is applied when studying convergence, i.e. estimating 



(^h m I? m )vP-n (*, Wo ,2/ 



Since Chb{r)t < 1 for h < h and t < t , ho, to sufficiently small, we get 

n/.N^f /mi./ ^ nris (?Chb(r)t) m (2Chb(r)t) m 

C(h)2C > (Chb(r)tYgt < C(h) - /n „, \ , \ : e h& < - , , { z e ^ 

v i \ > j at - \ J-i _ (2Ch b(r)t ) m ~ 1 - (2Ch b r )t ) m 



£=m 



where the ratio tends to zero for m — > oo. Moreover, we also get that C(h)f YlT=m(Chb{r)i) i < 

C ( h ) / 1 -(2ChoUy)t ) m ^ en< ^ s *° zer0 f° r 171 00 • Hence the series h m I® m con- 
verges uniformly. Hence we have shown that the solution u of the homogeneous 
system is given by a convergent series on any compact set if h is small. 

Remark The result is local in as much as u D is localized in D. But the Q in the 
above estimates in Theorem 14.21 is uniform. Since the manifold is compact this result 
is enough for getting global results. This concludes the studies of the heat kernel. 
Our result is related to the construction of a generalized heat kernel on a bundle 
over a compact manifold in |BGV91| concedrning the structure of the underlying 
space, however the diffusion is degenerate. 



5 A trace formula 

In this last section we use the uniform bounds derived above to derive a trace 
formula, i.e. an asymptotic expansion for the trace of the semigroup describing 
the evolution of the solution to equation ([]]). Let us explain the most important 
facts. Let H and H' be (separable, infinite dimensional) Hilbert spaces, and choose 
orthonormal bases (ej) i6 N an d (%)jeN m H and if', respectively. For any bounded 
operator A : H — > H' the quantity ||^4|| H 5 = Y^i j {-^- e h e j) i s independent of the 
choice of bases in H and H f , c.f. [Sim79||Roe98j. An operator if it is finite such 
that 1 1 A | \ HS < oo is called an Hilbert Schmidt operator and is called its 

Hilbert Schmidt norm. The Hilbert Schmidt norm is induced by a scalar product 
(A, B) = Ylii j (Aei,ej)(Bej,ei). A bounded operator is of trace class if there are 
Hilbert Schmidt operators A and B on H with T = AB. Its trace Tr(T) is defined 
to be the Hilbert Schmidt inner product (A*,B)hs- The trace Tr(T) depends in 
fact only on T i.e. it is independent of the choice of A and B. 
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In |BGV91| chapter 2 a covariant approach to a generalized Laplacian with corre- 
sponding heat kernel on bundles over a compact orientable Riemannian manifold M 
is developed. In particular, the scalar product between the space of compactly sup- 
ported sections of the bundle £ and the space of sections of £* <g> |A| is established, 
where |A| is the density bundle. The result given in Theorem 12.60 in [CFKS87] can 
be adopted. Since the generator C given by (Q]) considered in this paper is nonsym- 
metric we will, however, mainly exploit the existence of a kernel instead of working 
with the operators. Since the kernel is smoothing on the space of distributions ac- 
cording to |Hor68| and |01e71| and since it constitutes the transistion probability of 
a diffusion the operators exp[i£], t > 0, the corresponding Hilbert Schmidt norm is 
given by 



\e tc \ 



HS 



u{x 1 y 1 x 1 y)dxdy< I I u D {x ) y,x,y){l + 0{t))dydx 

T*M JM JTM 



< vol(M) / er^y (1 + 0(t))dy < oo 

JTM 

is finite, see e.g. c.f. [Sim79j. Exploiting the semigroup property of a transition 
density we see that the operators of the semi group are in fact of trace class 

Trexp[tC] = / u(t,x,y,x,y) dxdy = I u(—,x,y,x',y')u(—,x' ,y',x,y) dx'dy'dxdy 



■2' ' v 2' 



(exp[-£],exp[-£]}. 



Inserting the expansion we derived in Section 4 we find: 

Theorem 5.1. Under the hypothesis of Proposition l3A\ we have the following asymp- 
totic expansion with respect to a small parameter t. 

r 12? 

Trexp[t£] = / u(t,x,y,x,y)dxdy = - — -^-(volM+Rt 3 +a 4 t 4 +- ■ ■+a v t v +0(t v+x )), 

J T*M [2ix) n t ln 

where volM is the volume of M and the first non trivial coefficient R is the scalar 
curvature. All further coefficients dj, j = 4, . . . are invariant quantities of M which 
can be computed by the method indicated before. 
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